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SOME CONSIDERATIONS ON THE STRUCTURE OF 
QUADRUPED EXTREMITIES 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


An attempt is made to develop a general theory of the structure of . 
an extremity of a quadruped or biped, based on its mechanical function. 
The problem is essentially one of dynamics of a chain of linked levers. 
The equations of motions of such systems are practically insoluble ex- 
actly at the present time. A very crude approximation method is there- 
fore outlined. It is shown that, theoretically, the number of links in the 
chain forming the extremity may be determined by the optimum use of 
the total available force of the flexcrs and extensors. Discussion of the 
formula derived from this point of view shows, however, that actually 
this is not the determining factor. Other theoretical possibilities are dis- 
cussed. While no application can be made at present, yet in principle it 
is shown how mechanical consideration can throw light on the anatomical 
structure of the extremity. 


In a previous publication we derived a relation between the maxi- 
mum speed of locomotion of an animal and the length of its extremi- 
ties. We have found physical reasons for the necessity of the extremi- 
ties to be “folding” and to consist of at least two parts. Actually 
the structure of the quadruped extremity is more complex, although 
in most cases it may be considered within a rough approximation to 
consist of three parts: the hip, the ankle and the foot. In man the 
last is much shorter than either of the first two parts, yet in running 
the foot is still used as a lever, and not merely as a supporting sur- 
face for giving better contact with the ground. In some animals, as 
the elephant, the size of the foot is also rather small, compared to the 
size of the hip and ankle. In most other animals, however, especially 
in such as the horse, the dog, and the camel, the foot is as long as or 
even longer than that of the hip. Thus in the most general case the 
extremity may be roughly considered as consisting of three linked 
levers. Inasmuch as a system of two linked levers is sufficient to pro- 
vide for a “repelling” mechanism, one of the problems which we con- 
front is to find a mechanical reason for the general existence of three 
lever chains. 

Before we can hope to reach any conclusions that are applicable 
to reality, we must again make a purely theoretical study of different 
possible situations. This paper is the first step towards such a study. 
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A general dynamics of chains of levers was developed long ago by 
O. Fischer (1893) [cf. also Winkelmann and Grammel, 1927]. While 
the equations of motion of such a system are readily established, their 
integration even in simple cases presents almost insuperable difficul- 
ties, and recourse must be had to very drastic approximations. 

In attempting a dynamical theory of the action of an extremity 
as a mechanism which acts to repel the body from the ground, we 
must remember that strictly speaking we cannot consider an extrem- 
ity alone. We have to consider the whole organism as one system of 
linked levers. For such a system each equation of motion contains 
terms involving all other members of the system. Perhaps a possible 
simplification of the problem may be obtained by considering the over- 
simplified two-dimensional system represented on Figure 1. Here 


A 8 


ae 


c re) 
FIGURE 1 


the member AB plays roughly the part of the trunk, while AC and 
BD represent correspondingly the rear and front extremity, each con- 
sisting in general of n levers. Some fundamental properties of the 
corresponding real three-dimensional systems may probably be ap- 
proximately explained by such a schematized picture. Mathematically 
it has the great advantage of being a non-branching chain of levers, 
while actually any animal is a branching chain. 

However, even an approximate treatment of such a two-dimen- 
sional system represents still rather great difficulties. Therefore, we 
shall make here a further simplification, suggesting subsequently a 
possible method for handling the system represented in Figure 1. 

We shall schematize our picture further by considering only one 
extremity, consisting of m links, and by assuming that the mass M 
of the whole organism, minus that of the extremity, is concentrated 
near the upper end of the uppermost link. We shall number the links 
consecutively, beginning with the lowest, whose lowermost point is 
in contact with the ground (Figure 2). During the repelling action 
of the extremity this lowermost point P, remains fixed. 

We shall denote the mass of the ith link by m;, its length by J; . 
The position of each link is given by the angle ¢:, as indicated in 
Figure 2. 

We shall make further the following simplifications. We shall 
assume that each joint is actuated by muscles belonging to the near- 
est proximal link, that is, the ith joint is actuated by museles belong- 
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P, 


ing to the (¢+1)st link (Figure 3). Actually the moment of the mus- 
cular force with respect to the axis of the joint varies with the angle. 
which the two links make at the joint (O. Fischer, 1893). We shall, 
however, assume that this moment is independent of the angle and 
equal to the product of the constant muscular force F;,:., of the ex- 


FIGURE 3 


tensor muscle i + 1, by the constant shoulder 6, which we assume 
the same for all joints. All the torques b F;,i. act in such a way, that: 
they tend to rotate the links with ¢; > 2/2 clockwise, those with 
¢i < 2/2 counterclockwise. The latter direction of rotation will be 
considered as positive. Moreover we shall consider that at the point: 
P, of the nth link, in which the mass M is concentrated, a torque Fb 
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is applied, tending to move the mth link clockwise if m is odd, and 
counterclockwise if n is even. 

Following O. Fischer (1893) we shall introduce the following 
notations. We shall call the ith principal point H; the center of grav- 
ity of a system, which is obtained by concentrating the masses of all 
links, proximal with respect to 7, including the mass M, at the ith 
joint, the masses of all links distal to i at the (i—1)st joint, and the 
mass of the ith link at its center of gravity. We shall call such a sys- 
tem the ith reduced system. Denoting by S; the center of gravity of 
the ith link, we shall introduce the quantities c; and d;, which are 
defined in Figure 2. 

If, as we assume 


then with the exception of the nth joint, perhaps, the principal points 


of each link will lie very close to the proximal end of that link, and 
thus we have approximately 


Cj = l; . (2) 
Furthermore, we define the quantities 4; as follows: 4; is the 


radius of inertia of the 7th reduced system, with respect to the axis 
passing through the 7th joint. We also put 


M.D my=My« (3) 
Because of inequality (1), we have approximately 
Ma Mau: (4) 


Let @;, be an element in the ith row and kth column of the sym- 
metric matrix 


Pe ee omen, cane - Sand BET EES) 
eae ee Ae ® oS ea eS. 
BGs: Slate dae tie liars een rane ee (5) 


TO Ce a EFS Eager I: 
and let 
B,=—M, 9 ¢; cos ¢;. (6) 


Consider.a deformation of the system, in which the ith link is 
rotated by the angle d¢; , while all distal links remain stationary and 
all proximal are subject to a pure translation, so that none of the 
angles ¢;, & S iis changed. Then Bjd¢; represents the work done in 
such a displacement by the force of gravity. The work done by the 
torque b F;.,,; in such a displacement is 
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(—1) i+ b Fee dedi ; 
the work done by the torque b F;,is; is 
Cd Ltt 0 Me fae ey a (7) 


The work done in the displacement dd; by all other torques b Fy; is 
zero. 
Denote 


k=n d2 dk d¢ 2 
Ai=M, > aix[ecos (¢i — dx pace 
0B dxLeos (gi — gx) Saha 
Then the equations of motion of the system here considered are 
(Fischer, 1893) : 
Be Bye (1) OA haa 0 Fe) : (9) 


+ sin ($i — #1) ( 


Let 
i — Se ae (10) 
then, by a similar approximate argument, as used before (Rashevsky, 
1944), the time 4¢ during which the extremity exerts its repelling ac- 
tion is of the order of magnitude of 
At * L/v, (11) 


where v is the velocity of locomotion. 
Let, during At the angles ¢; vary on the average by the amounts 


die Then the average angular velocities are of the order of i / At = 


¢: v/l, and the average accelerations are of the order of ¢; v?/l?. 

In the final extended state of the extremity the values ¢; — 9x 
all become zero; hence, the cosine terms become equal to one, and the 
sine terms vanish. If we consider the case, when initially the values 
of 6: — ¢x were not too large, we may very roughly consider during 
the process of extension all cos(¢: — ¢x) terms as well as all 
sin(¢; — ¢.) terms as constant, with values between zero and one. 
We may do the same thing for cos ¢; in the expression for B; . 

Putting 

cos (¢i—¢x) =an; sin(¢i—¢dx) =Bix; cosdi=—yi; (12) 
and averaging the equation (9) by substituting 


‘pias v/l; oa Bt wrt, (13) 
we obtain the following system 
M, v? ot me 
: Dd dix (ain oi + Bix 624) =— Mog Ci yi + (aed) CD ae 5 
k 


(14) 
a b FG as) 
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If all forces Fiji are prescribed and the geometrical structure, 
that is, the J;’s and the c;’s are known, equation (14) gives n equa- 


tions, quadratic in ¢;, from which the values ¢; can be determined. 
Thus, given the initial position of the extremity, we may roughly de- 
termine final position after the time 4t = v/1. 

It may be more interesting, perhaps, to reverse the problem and, 


considering all the $i as given, to determine the geometrical struc- 
ture which minimizes the sum of all the muscular forces that are nec- 


essary to produce a given displacement ¢; . 
For i= 1, equation (9) gives 


b Fy = (—1)"? (A, + B). (15) 
Substituting that into equation (9) (with i= 2) we find 
b Fog =: (—1) *° (A, + Ao + B, + B.). (16) 


By substituting this into the equation (9) with 7 = 8, we obtain an 
expression for b F;,, and thus by successive steps we obtain in gen- 
eral: 


b Fein = (—1) 4 (A, + Aa +--+ As +B, + Bp +--+ + Bs), (17) 
and finally 
b Hiss. (a1) (Aj h:*shAs # By tee th Ba)s (18) 


The geometric sum of all the torques is readily obtained from 
equations (17) and (18) and is found to be, 


20 Fi 40 = 3 SS" LC Ee) (A Ba) etl 


It:is; however, the sum of the absolute values of all the torques b Fi,is: 
that is of more interest, since it is the latter that determines the size | 
of the necessary muscles. It is equal to 


bF,=> |b Fy in| =n(A, + By) ss 
(n—1) (A, + B,) 5 (A, + Bak s 


With the assumptions which we made and which lead to the approxi- 
mate equality (2), the A,’s are quadratic forms in the J,: 


Dfekh, (21) 


(20) 


with 
fin = (ain bi + Bix 6s). (22) 


The B;’s are linear in the 1;’s. Hence the right side of equation (20) 
iS a general quadratic polynomial in the J;’s with linear terms. We 
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may ask for values of |; which make that polynomial a minimum with 
the condition (10). This for a given n, would give the relative sizes 
of the different parts of the extremity. The problem reduces essen- 
tially to a system of n linear equations in the l,’s, but its general 
treatment is rather messy for arbitrary n and shall not be given here. 

Instead we shall assume that all l,’s are for some reason equal 
or approximately equal, so that 


L=-. (23) 


Each term of any A; varies, therefore, as 1/n?, and since there are n 
terms in each A; , therefore, each A; varies as 1/n. The B;’s vary as 
1/n. Denoting by h; the coefficients of proportionality, so that A; + 
B; =h;/n, we have from equation (20) 


1 
Oia spi Se tiind sae 1 high Sc: he) - (24) 


If h* is the average of all h;’s we have approximately 
mgs (1 +---n) iD) E 
n 2 
Hence F', increases indefinitely with increasing n , and from the point 
of view of the minimum force, an extremity consisting of only two 
parts is the best. 

Let us now consider a different situation, namely, when the ex- 
tremity is lifted by the flexor muscle to be thrust forward during the 
‘flight phase” of running. Now we may consider as a first approxi- 
mation the position of the center of gravity of the whole organism as 
prescribed, and study the movement of the extremity with respect to 
the center of gravity. In this case we may consider the extremity as 
having a fixed point P; at its upper end. Since now the movements of 
the extremity do not involve the displacement of the mass of the 
whole organism, we do not need to concentrate the mass M at the 
upper end and only consider the masses m; of the individual links. 
Considering those approximately as homogeneous cylindroids, we may 
put 


(25) 


Mix L; é (26) 


We now shall number the links from above downward, the first being 
thus the uppermost. Further we must remember that the distances 
l; are measured from the principal point of the ith link to the joint 
which lies in the direction of the fixed point (now P,) of the chain 
(Fischer, 1893). Therefore now the c¢;’s will not be even approxi- 
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mately equal to the corresponding i,’s. If all l,’s are equal and if 
equation (26) holds, then c; decreases with increasing a. For any 
l,’s, ¢; is proportional to l;, the coefficient of proportionality being 
a function of the index 7 and of all the others , k = 71. The equations 
of motion of the system are still given by equation (9), but in the 
expression for A; and B; we must substitute for M, the value 


Mo => Mi. (27) 
Because of the proportionality of the c;’s to the l;’s, the expression 
for A; reduces again to quadratic form in 1;l, of the same type as be- 
fore. The only difference now is that bF is applied to the 1st link and 
not to the nth. In the same manner as before we now find: 


b.F p= Ay + Buch 2 (Asse Bsa st As aa): (28) 


and we see, that again for equal J;’s , fF, increases with n. 

There is, however, one case in which the situation changes. It is 
clear that if most of the weight of the extremity is located in the up- 
per part, it will take less work to lift the extremity during the flight 
phase. Such a distribution of masses is actually the case. Let us, 
therefore, consider the following case. Let all the links except the 
most proximal be geometrically similar, so that their mass would vary 
as 1;3. The thickness of the link may then be determined by con- 
siderations of bending strength of a rod of given length. If we 
now increase n, keeping m, constant, we shall have a shift of the 
center of gravity of the whole extremity in the proximal direction. 
The upper link will be much thicker than the other. With increasing 
nm, therefore, all H;’s will shift in the proximal direction, and there- 
fore the ¢;’s will not remain proportional to the /;’s, but decrease 
more rapidly. For equal /;’s every term in A; will now decrease more 
rapidly than 1/n?, and the whole A; will decrease more rapidly than 
1/n. If, as a very rough approximation, we put now ¢; « 1;?, then 
A; ~ 1/n’. The B;’s also vary as 1/n?. The expression for F’, now 
becomes approximately of the form 

h’ h' 


n(n + 1) 

h’ ak Sacteegee 

Qn? 2 Y 2n’ (23) 
where h’ is the average value of the coefficients of 1/n in the differ- 
ent A; + B;. Thus in this case F’, decreases with increasing n. 

For the repelling action of the extremity, when we again have 
equation (2), the old result holds, namely that F, increases with n. 
Hence, denoting by C a constant 


bF',4= 


, 


h 
Woh tarmteks pl Se Phere Sd (30) 
n 
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Hence bF’; has a minimum for 


ae 5 
nr— =a 
‘i 2h* au 


The value of h* is approximately 


h* © M,[v?(aix di + Bix 6:2) + gyll, (32) 
while 

hi © Mol[v? (ain di + Bix $2) + gyl). (33) 
If for the repelling action of the leg and for its lifting during the 


flight phase the average values aid; + Bindi? are approximately the 
same, then 


n= | hes 34) 

VM 
With m./M, © 107 and 1 © 10? em, n is about 3. But for shorter ex- 
tremities it would be less than one, hence only 2 links would exist. 
Thus while these considerations lead us to the determination of the 
optimum number of links of an extremity, they have only a theoreti- 
cal value, for in reality, obviously, that number is determined in some 
other way. 

It must be remembered that unless we make some rather artifi- 
cial assumptions about the attachments of the muscles, we should con- 
sider in equation (29) b itself as variable, namely put b = kl; = kl/n, 
k being a constant. The forces F', and therefore F’, also is likely to 
vary as b? = k°I?/n? , being of the form fb?, where f is essentially the 
same constant as used in our previous paper (Rashevsky, 1944, p. 28). 
The constant f is determined by the physicochemical considerations of 
the properties of the muscles and is probably one of the determining 
parameters of the species. Putting in equation (30) b = kl/n; 
F’, = fk’I?/n?, we obtain a fourth degree equation for the determina- 
tion of n. If this equation has more than one real positive root, the 
smallest must be chosen by the principle of maximum simplicity (Ra- 
shevsky, 1943). It is worth while to investigate, what values for n 
are found in this way. 

The reason for having usually three principal members in quad- 
ruped extremities may be also an almost kinematic one. For a re- 
pelling action two members are enough; however, if there are only 
- two members (Fig. 4a), then the following happens during the “un- 
folding” of the extremity. While the second member rotates clock- 
wise by an angle dd», the first member rotates counterclockwise by an 
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FIGURE 4 


angle d},. Both rotations tend to lift the body, but while the rota- 
tion of 2 moves the body to the right (forward), the corresponding 
rotation of 1 moves it to the left (backwards). As we have seen (Ra- 
shevsky, 1944) in fast running, the body must acquire a velocity, 
which forms a rather small angle with the horizontal, due to the re- 
pelling action of the extremity. Hence the horizontal component of 
the velocity imparted to the body must be much larger than the verti- 
cal component. This will not be the case for an extremity consisting 
of only two members. With three members (Figure 4b) the situation 
is different, for now two members rotate clockwise, and one counter- 
clockwise, increasing the horizontal component. With four members 
we again have the same situation as with two. 

For walking, which does not require much folding of the extrem- 
ity, and where the extremity rotates in the hip almost as a whole, 
except for a slight bending in the knee to reduce the lift of the center 
of gravity, two members may be sufficient. Elephants with very short 
metatarsus do not run, but walk rapidly in a rather rigid fashion 
(Howell, 1944). 

The mathematical treatment of this aspect of the problem would 
be as follows: Equation (14) should be solved for given Fi’s and 
prescribed 1;’s for the n average displacements ¢;, which in general 
will be of different signs. If the ¢;’s , as well as the initial values of 
¢i 8, are known then the horizontal displacement of the uppermost 
point of the extremity is also known. Thus we may calculate whether 
it is larger for n = 8, than for n = 2, and, if so, how much. 

We may consider the ratio of the total kinetic energy T to the 
kinetic energy T, of the rotation with respect to the center of gray- 
ity, and investigate for what geometrical structure of the extremity 
that ratio has a maximum. If we consider that n = 3 is the optimum 
from the above kinematic considerations, then we may determine I, , 
l, and 1; for a given 1 from considerations of the minimum of T,. The 
expression for 7’. is (Fischer, 1893, p. 32): 


T= 5 mo Pha? gh? + Bea* a? + Ks? ga + 2dh 2 COS (1 — $a) $s de 


~ (35) 
+ 2d, C3 COS (1 — $3) hi bs + dy C3 COS (bd2— $3) bo dol; 
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and for T, = T — T, (Fischer, 1898, p. 35): 


1 ’ é : ae 
gan 5 [¢,? pr? + Oy? bo? + Cy? bs? + 2C, Co COS(d1— $2) $1 $2 
(36) 
+ 2¢, Cs COS (¢1 — d3) oi $3 + 2Ce C3 COS (bo — $s) ba dol. 


The quantity k; is the radius of inertia of the ith reduced system with 


respect to H;. If the values of 4; , correspondingly of $iv/1, are pre- 
scribed, then the problem is to minimize a quadratic form in J;’s, re- 
spectively in the c;’s. This problem reduces to a set of the linear equa- 


tions. We may investigate whether for all ‘hi equal we have all J; or 
all ¢; equal, while a larger ry corresponds to a shorter 1; or c;. If 


such were the case, it may account for a shorter metatarsal part in 
plantigrade animals, since that part would have to be moved through 
a larger angle during extention. 

In all cases the solutions will depend on v and / as parameters. 

The author is indebted to Dr. A. S. Householder for a critical 
reading of the manuscript of this paper. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Fund of the University of Chicago. 
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TREE TRUNKS AND BRANCHES AS OPTIMUM MECHANICAL 
SUPPORTS OF THE CROWN: II. THE BRANCHES 


MARTINUS H. M. ESSER 
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Using assumptions made in a previous paper, a theory of the shapes 
of primary branches is developed. Two cases are studied: a primary 
branch which has a portion denuded of secondary branches, and a pri- 
mary branch with a continuous load of secondary branches. The first 
case leads to hyperelliptic integrals, the second to polynomials of 
second degree. : 


Denuded stem. We consider a branch as described in the intro- 
duction to a previous paper (Esser, 1946; hereafter referred to as 
Part I). The same notation will be used. The shape of its stem is 
determined by equation (4) of Part I. The first two terms of the 
right side of this equation have a ratio of the order of 

cos 6 T(z) 
sin 6(¢ — 2) ; 
When @ is not small, this ratio is small and we can therefore disre- 
gard the first term of the right side of equation (4) of Part [. Simi- 
larly, 
G cos 6r (2) 
ao 
can be neglected in comparison to 
4G 
— (Y—zsin@). 


ao 


Finally, we neglect also the action of wind. Equation (4) of Part I 
thus reduces to 


iné " 4G 
rz) =APEES fH eyee(endt += (¥—zsing). (1) 


Equation (1) determines 7 as a function of 2. 
We introduce the following notations: 
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i 4Gsin 6 iY. 
Annee Bae eee 
o Io sin 6 
Then equation (1) becomes 


yi =A [ C—avr(d +BX—2). (3) 


Two successive differentiations give 


d “1 
Via | peOat Ss, 
dz z 
a (4) 
4 Ay/(z). 
dz? 


We consider y’ = dy/dz as a function of y, and obtain from the pre- 
ceding equation 
d? dy d d d 
Os Se ES pe eas 
dz? dz dy dz dy 
d(y')?=2Ay*/*dy . 


Integrating and denoting by a the constant of integration 
6 
y? is Ay®/? + a. (5) 


From physical considerations it follows that the function 7r(z) is de- 
creasing; therefore y’ is negative. The preceding equation gives then 


6 
—A 5/3 + @ 
ea 


Integrating and denoting by b a constant of integration: 


= a2. 


} d 
»- f=, (6) 


" 6 
#o,[8a 05/3 + @ 
5 


We now determine the two constants of integration. Denote by p the 
radius of the cross-section where z=1. We have from equation (1): 


314G 
p=r(l) =, /— (Y—lsin8), 
ao 
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and from equation (4) we have 


y (1) =—B D) 
from equation (5) we have 
6 6 
a = B? —— Ay (l) = B?—- Ap’; (7) 
5 5 
and we find b = 1 by making z = / in equation (6). 
Thus we find 
(rae ‘ gsds 
z=l— — = 1 — | ——_., 
6 
eace e+e A (65/3 — p>) + B? p JA + B? 
Hence we obtain, as a final result, the following relation between r 
and z 
s? ds ne 
Reece 5 
wi A (s§ — - ee 
9 . 
where 


: 4q@ sin 0@ 4Gsin @ 
—lsiné), A=———_, B=——_.. 


oc ao 


Weight of branch. We will now find a simple expression for the 
weight M; of the branch. This expression will be used in a subsequent 
study. 

We have 


1 
Ms=o2 | re(t)dt + G. (9) 
A combination of equations (4) and (5) gives 
4v2(0) +a—B 
ip a praee PL Tera ee 3 
get Aa) OL A = x 


From equation (7) and from the last equation (2), we find 


6 BAe. 
froa= ls [rs(0) — eo 
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Substituting this value into equation (9) and using equation (2), we 
obtain the final result 


Maa | ee [r>(0) —7(1)] + G?. (10) 


10 sin @ 


We shall now study the particular case where G = 0. This case has 
already been commented on at the end of Part I. Equation (8) re- 
duces to 


(11) 


Substitution of this value into either equation (9) or (10) gives 
* 42 w' sin? 6 : 
Mere. 


The interest of this formula is to show that, in our approximation, 
the weight of a branch varies as the fifth power of its length. 

Continuously loaded stem. Instead of considering a primary 
branch as consisting of a stem terminated by a crown, we now con- 
sider a primary branch to consist of a stem loaded with a continuous 
distribution of secondary branches. Let p(z) be the load per unit 
length of stem in the neighborhood of z. Then the weight of the por- 
tion of primary branch (stem plus its load of secondary branches) 
between the cross-sections z and z + Az is 


Mo {fess v2 (t)dt + {Res p(t)ac. 


(12) 


The integral equation determining r(z) and corresponding to 
equation (1) of the preceding problem is now 


4 sin 6 


o 


re (2) = 


iP c— 2)| r2(¢) + a4 at. (13) 


Two differentiations give: 


d’rs(z) 4siné p(z) 
E r*(2) + : 


dz? n 
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We will solve equation (13) in one particular case, namely the case 
where the following four assumptions hold. 


I. The length of the secondary branches decreases to zero when’ 
we approach the tip of the primary branch. We shall consider a linear 
decrease, the length 1,(z) of the secondary branch attached at the 
point z of the primary branch being given by l,.(z) « (l— 2). 


II. The number of secondary branches increases as we approach 
the tip of the primary branch. We shall consider 


1 
n(z) Sain 7" 
i—¢ 
where n(z) is the number of secondary branches per unit length of 
primary stem in the neighborhood of z. 


III. The weight of the secondary branches is proportional to 
the fifth power of their length, as in equation (12), 


M2(z) « 1,5(z). 


IV. The secondary branches can be grouped into pairs of branch- 
es symmetric relatively to their common point of attachment to the 
primary stem. 

Assumption IV allows us to replace each pair of branches by a 
weight located at their point of attachment. We obtain from the first 
three assumptions, 


p(z) =2(z)M2(z) « (l—2z)*4, 
p(z) =A(l—2z)!, 
A being some coefficient of proportionality. 
Equation (13) now becomes 


4 sin 6 


if (¢—2)| or*(©) += (i= 2)+ Jae. 


(Zz) = 
By substitution, we see that this equation has the solution 
r(z) =a(l—z)?, (14) 


where a is the solution of 


2 sin 6 A 
a= “ite ap hy 
105: PL 


oO 


This equation has one and only one root a. This root is positive. (This 
result can be seen by making a graph of both members of the equa- 


tion.) 
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Equation (14) determines the shape of the primary branch. 
This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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The calculation of the size of the “sensitive volume” or “control 
center” in biological effects of radiations is discussed from the viewpoint 
of the probabilistic theory of these phenomena based on the concept of 
random “effective events.” On the bases of that theory, the resistivity of 
a microorganism to radiation is defined as its “mean life” under a radia- 
tion of one roentgen per minute. This mean is calculated for processes 
with and without recovery. The case of variable sensitivity, as it oc- 
curs for instance during mitosis, is discussed in detail. Methods are 
given to calculate this variability from survival curves or similar experi- 
mental data. The theory is applied to experiments of A. Zuppinger on 
irradiation of Ascaris eggs with X-rays. 


1. The physicochemical background of the theory. The mech- 
anism of biological effects of radiations is complex. In the field of 
X-rays it is agreed that one of the primary effects of radiation is the 
ionization which is produced either directly by the radiation quanta 
or indirectly by the electrons ejected from the atoms. This ionization 
is in general only a first step in a chain of phenomena leading to the 
observed effect and consisting in general of some stable chemical 
change (See e. g. Clark, 1940, chapters iv, viii, x, xi; Clark, 1944; 
Failla, 1936). Such a change involves, according to the present con- 
ceptions, some vital control centers of the microorganism. A center 
of this kind may consist of a very small aggregate of molecules, or 
of a single molecule and perhaps even of an aggregate of atoms in- 
side a molecule. As P. S. Henshaw (1944) points out the functional 
activity of a cell involves a relatively large reserve of material, which, 
if damaged by radiation, may be renewed more easily than those sub- 
stances on which the control activity of the cell is based. These latter 
substances are in general present in much smaller amounts, and it is 
more difficult to restore them to their original condition once they 
were inactivated; besides this, their action seems to be dependent to 
a high degree on their concentration, so that even changes involving 
small amounts may produce substantial effects. These and other theo- 
retical reasons based on a mathematical analysis of experimental] data 
(Opatowski, 1945b, pp. 162-163) suggest that biological effects of ra- 
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diation may have a primary relation to the changes occurring in an 
extremely small part of the microorganism. 

The theory of the biological effects of radiations has been limited 
up to now to a probabilistic approach to these primary processes con- 
sisting of ionizations and dependent directly on the random absorp- 
tions of radiation quanta by the control centers or sensitive volumes 
of the organism. However it must be kept in mind that what actually 
is observed is an effect of the tonization. This effect may be reached 
in three different ways: (i) through some action which is developed 
practically instantaneously; (ii) through a chemical or biochemical 
reaction, and therefore through a process of a deterministic charac- 
ter requiring a certain time after the absorption of the radiation 
quanta occurred; (iii) through a combined action of both mechanisms. 
The latter is the most general type of process involving a superposi- 
tion of random and deterministic phenomena. A mathematical analy- 
sis of some experimental data suggests that this general mechanism 
may actually occur. The second process representing the so-called de- 
layed action is substantially identical with the first, since it involves 
only a translation of the origin of the time axis. The present theory, 
as well as those previously published, concern the first two mecha- 
nisms. These theories are based on the concept of an effective event 
occurring in a sensitive volume or control center of the microorgan- 
ism. Such an event may be the absorption of a photon or an ioniza- 
tion; it may cause the formation of a new substance inside the sen- 
sitive volume and elsewhere, or may change other physical charac- 
teristics of the cell leading to such phenomena as cell swelling, 
variations in concentrations and in rates of metabolism, etc. In 
the mathematical theory the occurrence of an effective event in the 
sensitive volume is considered as a transition of that volume to a new 
state, The observed effect is reached after a certain number of transi- 
tions between the successive states has been carried out, so as to reach 
a state n which represents the observed effect (cf. Opatowski, 1942; 
1945b; 1946). The probabilistic model used (Markoff chains) has 
been defined in a previous paper (Opatowski, 1946). It is the same 
model which has been introduced by E. Rutherford into the field of 
radioactive transformations, except that we consider here also the 
possibility of reverse transformations representing the phenomenon 
of recovery. The recovery may be due to reversibility of reactions, 
to occurrence of new reaction chains or to ionic recombinations. The 
states are numbered progressively 0,1,---,. During an infinitesi- 
mal interval of time dt transitions between two consecutive states 
are considered as the only possible, which is a limitation of the 
theory. The probability of a transition (i—1 — i) within an interval 
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of time (¢---t+dt), under the hypothesis that at the time t the state 
t—1 existed, is assumed to be kidt , where k; is not necessarily a con- 
stant. A similar probability for the transition (i+1 — i) is put equal 
to gidt. Under these assumptions the probability Y,(t) that the 
microorganism reaches the state n at a time ¢ of irradiation of a 
sample of microorganisms has been calculated (Opatowski, 1942, 
1945, 1946). The function Y,(t) is the global transition probability, 
inasmuch as it represests the probability of transition from the initial 
state to the state which is under observation. (A similar probabilistic 
model has been suggested recentiy also for the action of viri and anti- 
bodies [see Jordan (1941) ; v. Schelling (1942) ; Koyenuma (1942) ]. 

From a formal viewpoint the theory may be interpreted as rep- 
resenting a kinetics of a chain of reactions, which in general may be 
reversible. However, the interpretation based on the concept of ran- 
dom effective events gives a possibility of calculation of the size of 
the sensitive volume. For this purpose the effective events are inter- 
preted as absorption of photons or formation of ions. For instance, 
in the case of X-rays, the calculation is based on the assumption that 
the energy absorbed by the irradiated sample is measured in roent- 
gens. This measurement gives that fraction of the total radiation 
energy E which would be absorbed in ionizing 1 cm’ of air. Conse- 
quently the number of ions produced in 1: cm* of air by the energy 
E may be calculated from the definition of the roentgen (2. 1 < 10° 
ion pairs per roentgen in 1 cm of air (cf. Clark, 1940, p. 192). The. 
number of ions produced is proportional to the absorption coefficient; 
this enables us to calculate the number n, of ions produced by FE in 
1 cm? of irradiated matter, as soon as the chemical composition of the 
latter is known. A comparison of the experimental results with the 
mathematical expression of the probability function Y,(t) gives the 
number of effective events n produced in the sensitive volume by the 
energy E. Let » be the number of ions in the sensitive volume corre- 
sponding to one effective event. Then »m/n, is the size of the sensi- 
tive volume. At the present stage of our knowledge of the action of 
X-rays, only a range of values may be assigned to »; consequently 
what the above calculation gives is only a lower and an upper limit 
for the size of the sensitive volume. This fact does not seem to have 
been noticed by some authors. If formation of each ion, indifferently. 
of its sign, is an effective event, then » = 1. If ions of a certain sign - 
count only for the biological effect, then y = 2. If collisions of the 
primary photons with the sensitive volume are the only effective 
events, then v may be very large within wide limits, since an electron 
of 50,000 eV may produce 1,000 ions, whereas one of 2,000,000 eV 
may produce more than 100,000 ions. Comparative experiments with 


104 EFFECTS OF RADIATIONS 


different wave lengths could perhaps give some information as to 
the nature of the quantity ». 

It should be, however, kept in mind that the above method 
of calculation of the sensitive volume is based on several approxima- 
tions: First of all, the irradiated microorganisms are physically non- 
homogeneous, even if they form a biologically homogeneous aggre- 
gate. The X-ray beam, even if monochromatic at the source, may 
become polychromatic inside the living matter. The formation of ions 
along the electronic track is not uniform, being more dense at its 
end, because of lower electron velocities. Elements of higher atomic 
number are loci of intense ionization, since the absorption coefficient 
increases rapidly with the atomic number. The fact that these ele- 
ments form only a very small fraction of the irradiated matter does 
not preclude their possible importance for the final effect. In fact, it is 
well known how essential small amounts of impurities are for cer- 
tain photochemical phenomena (cf. e. g. Hill, 1946; Noddack, 1926, 
p. 605). A theory which would take into account all these facts would 
be complicated, although some attempts in this direction have already 
been made, for instance, the theory of R. Glocker, (1932) which does 
not neglect the possibility tnat an ionization of the sensitive volume 
may be due to an absorption of a radiation quantum outside of that 
volume. In order to consider the possibility of occurrence of a pri- 
mary effective event outside of the sensitive volume, some authors 
suggested recently the concept of a physically sensitive volume to be 
distinguished from the actual biological sensitive volume (Méglich, 
Rompe and Timoféeff—Ressovsky, 1942). However, as far as the 
ionization is concerned, this possibility was already taken into ac- 
count in a quantitative manner in the theory of R. Glocker. 

The physicochemical background of the theory of biological ef- 
fects of radiations has been outlined here in a very concise manner: 
for more detail the reader is referred to the literature already quoted 
as well as to the following publications in the field of X-rays and their 
applications: Compton—Allison, 1935, chapter vii and pp. 5-18; May- 
neord, 1938, 1940; Kirchner, 1980; Gray, 1936; Victoreen, 1944; 
Behnken, 1926; Lauritsen, 1933; Webster, 1933 ; Condon—Terrill, 1927; 
Glocker, Hayer and Jiingling, 1929; Scott, 1987: Braun and Holthu- 
sen, 1929; Grebe, 1926; Rump, 1927, p. 290. 


2. The object of the paper. In certain cases the seusitivity of 
the organism to the radiation is variable. This occurs particularly 
during mitosis, and must be taken into account when the duration of 
the irradiation is of the same order of magnitude as the life cycle of 
the microorganism. An increase of the sensitivity to radiation may 
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appear in two different forms: (i) The number of effective events 
necessary to obtain a certain effect may decrease; (ii) the probability 
of occurrence of an effective event may increase for instance because 
of an increase of the sensitive volume. In the first case the number 1 
changes, in the second case the quantities k;, g; which express the 
probabilities of transition from one state to another change. The first 
case does not present any mathematical problem: it just means that at 
certain moments the number n changes by 1, that is, that the global 
transition probability Y,,(¢) has a finite jump at that moment. The 
second case means that k; and g; are functions of the time. We con- 
sider this case in detail, mainly under the assumption that no recov- 
ery occurs (g; = 0). The theory has as its objective the determina- 
tion of the functions k;(t) from experimental values of Y,(¢). Such 
determination is carried out explicitly under the assumption that all 
the k;’s are equal to each other; this involves either the approximation 
of taking for the k;’s a suitable average value k, or the assumption 
that an inactivation of only a small part of the sensitive volume is 
sufficient to produce the observed effect. The actual calculation of the 
function k(t) involves a first trial computation on the basis of an 
arbitrarily chosen value of 7, which is then improved to the correct 
value. The theory indicates the type of experiments which are neces- 
sary for the calculations. The computations, which involve the use 
of the incomplete gamma function, are applied to a set of experiments 
on irradiation of Ascaris eggs with X-rays. 

A new mathematical formulation of the resistance of a micro- 
organism to irradiation is introduced. It is based on the concept of 
average duration of global transition, i. e., of the average time neces- 
sary for the transition (0 > n). Such time L represents the average 
life of the microorganism in the case in which the state 7 represents 
its death. In the case in which the intensity of radiation is 1 roent- 
gen/min, L is measured in roentgen and is defined as the resistivity 
of the microorganism. General expressions for the calculation of L 
are given and explicit formulae are obtained in the cases in which the 
k’s and g;’s are independent of 7. 


3. The average duration of global transition (the mean life). 
We consider here only the cases in which the nth state is final, that . 


is, Kn =O. 


(i) We discuss first three types of processes in which the tran- 
sition intensities k, and g, are independent of the timet. 

(a) Consider the case in which no recovery occurs, that is, the 
intensities g, of all reverse transitions are zero. Then Y,(t) is a 
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monotonically increasing function from 0 to 1 with Y,(«) =1 (Opa- 
towski, 1942, pp. 84, 87). The mean life of a microorganism between 
the states 0 and n is defined as the mathematical expectation of f¢: 


ee [ tav.= [evecat. (1) 


Deterministically this definition is intuitive, because dY, is the frac- 
tion of the total number of microorganisms which reach the state ” 
between t and t + dt, i. e., which “Jive” a time ¢. (The concept of 
average life or of average duration of transition between the states 0 
and n is analogous to the concept of mean life of an atom in a quan- 
tum state in theoretical physics, (e. g. Richtmyer and Kennard, 1942, 
pp. 304-305). From a formula concerning the moments of the func- 
tion Y’,(t) (Opatowski, 1942, p. 87) we have 


L=31/kk=n/k, (2) 


i-1 


where k is the harmonic average of all the k;’s. 

(6) Consider now the case in which the reverse transitions occur 
between intermediate states only, that is, a process in which all the k,’s 
and g,’s are different from zero, except ky: = Gn1 = 0. It is known 
(Opatowski, 1945b, pp. 170-174) that such a process is equivalent to 
a process consisting of direct transitions only between the same num- 


ber of states with suitable transition intensities k,,--- , k,, n being 


a final state also in this new process. The quantities —k; are the n 
roots of the equation in s: 


[@,,¢||n=0, (3) 
where 


Arr =8 + Kina tr Gris Ara = k,; Oy ra = Gr; 
a,,.=90 for |r—e|>1; (4) 


0Sr7r,cSn—1. 
Consequently L = S‘=" 1/k;. But (Opatowski, 1945b equation 28): 
iE je oa ky = &, ke ++ Ky} 


L=ay4/ (k, Ka+++ Kn), (5) 


where @,. = © k, k, «++ is the elementary symmetric function of de- 
gree n—1 of the k;’s, that is, a,_, is the coefficient of s in equation (3). 
Consequently (cf. Opatowski, 1945b, equations 15 and 16) a, is the 


therefore 
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sum of all the minors of order n—1 of the determinant (3) (4) cal- 
culated fors=0. 

Consider as an example the case in which the k;’s and g;’s are 
independent of 7, that is: 


Then, the following explicit form of the equation (3) (4) may be used 
(Opatowski, 1945a, equations (17), (6): 


E,—gke—0, (6) 
where £,, equals (kg)”/? multiplied by a Tchebychey polynomial in 


o=(s+k+g)/(2Vkg), 


that is, 
E,= (kg)"?[(¢ + Vo? —1)™ — 
eee es Os cian. (7) 
ie hen Lee aot 1). 
The coefficient of s in the expansion of EF, is easily calculated as 
(dE,/dS) so = kh” (1—7)?[(n + 1) (1 +9) — 
(8) 


ee ot RG Baad hal Od cad 


where r = g/k. Consequently, dn... the coefficient of s in equation 
(6) is: 
he 1) fa (0 1 ee] 


and from equation (5) we obtain the average life: 
L=[n— (n+1) r+7r)/[k(1—1)?] = 
(n/k) [1+ (1—n")r+.---]. 


If r = 0 this expression reduces to (2). If there is recovery, the aver- 
age life increases, in a first approximation, proportionally to the ratio 
r of the intensities of direct and reverse transitions. 

(c) Consider now a process (p) in which the nth state is final 
and a reverse transition from the nth state is possible, that is, kn=0, 
- Ona # 0. As usual, let Y,(t) be the global transition probability 
for this process. Consider the equation 


||@r,c||mu =O, 0Sresn, (10) 


where the determinant ||@,,-|| is defined in the same way as in equa- 
tion (3), except that its order is now one unit higher. We know (Opa- 
towski, 1945b, pp. 174-175) that the relation (10), considered as an 


(9) 
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equation in s has one root equal to zero, and all the remaining roots 
negative numbers —&, —kz , +++ , kn. Let Yn(*) be the global tran- 
sition probability for a process (p) consisting of direct transitions 


only between the contiguous states 0,1,---,n with the intensities 
respectively k,, Fs, ++:, kn. Then (Opatowski, 1945b, pp. 174-175): 
Yn (t) = hy k++ hin Yn(t) / (Ba hee -++ Kin). (11) 


Since Y,,(«) =1, we have now (Opatowski, 1945b, p. 776) : 
Yn(00) = hy hy +++ Ken/ (Her K++ Hen) <1, 
and equation (11) may be written 
Y,(t) =Yn(0)Yn(t). 


To represent the average life, the formula (1) must be modified in 
the following manner: 


Yn (oo) 1 — 
L= p tdYn/Yn (co) = { tdY, . 
0 0 
We see in this way that the average life in the processes (p) and (p) 


1s the same. Therefore: 
L= Sif. (12) 


tl 
If we divide the determinant (10) by s, we obtain a polynomial: 
R= by + Ana S +--+ 8", 


and by known properties of the elementary symmetric functions of 
zeros of polynomials we have from equation (12) 


L= Qy-4/b, « 


As an example consider again the case in which the k;’s and gis 
are independent of 7, that is, 


k=k(1SiSn); 9=g9(0 StSn-1). 


Then (Opatowski, 1945a, formula 5): R = E,. Consequently d,-1 


yee the expression (8) and b, is the value of the expression (7) 
or s=0: 


bp, = "(1 — r™) /(1—1). 
We obtain in this way: 
L=[(m + 1) (1m) (1 + 9) — (1-7) 2 (14 1) 1/[k—)] 
= (n/k) [1 + (1—2n) r+..-], 
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Also here the recovery increases the average life in a similar manner 
as in the process (b). 


(ii) We consider now a case of variable sensitivity without re- 
covery, that is, without reverse transitions, so that all gj’s are zero. 
We assume that the variability of sensitivity is reflected in a depen- 
dence of the k,’s on the time in the form 


ki =K; f(t), (13) 


where the K;’s are independent of ¢. If we introduce a new time 
variable 


T= | f(tyat, 


we obtain the probability Y, by taking its expression valid for the 
case (ia) and changing in it formally ¢ into JT (Opatowski, 1945b, p. 
165; 1942, p. 84). However the calculation of the average life must be 
made with reference to t as the time variable. This calculation can 
be carried out explicitly when ¢ is a polynomial in T: 


ba ad 


and if T — + © whenever t ~ + ©. Under these assumptions, and 
choosing t = 0 as the lower limit in the integral which defines 7, we 
can write: 


1 ao 
Le { || tY’, (1) dT = 
0 0 
(14) 
Seo f T™ Y’,(T) dT =Sm M! Cm bn (Ket, +++, Ky), 


where h,,(X1,°-: tn) is the complete homogeneous symmetric function 
of the x;’s, that is, the sum of all the possible products of the 2;’s of 
degree m , each of the x;’s being taken any number of times up to m. 
(These functions are called also homogeneous product sums or alef 
functions of Wronski.) The passage from the fourth to the fifth 
member of equation (14) involves the use of a known relation con- 
cerning the moments of the function Y’,, (Opatowski, 1942, p. 87). 


4. Determination of the change of sensitivity of the microorgan- 
ism from experimental data. Since processes with recovery are re- 
ducible to those without recovery (Opatowski, 1945b, pp. 170-177) 
we may assume the intensities of reverse transitions g; to be all zero. 
We consider the intensities of direct transitions k;(t) to be functions 
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of the type (13). As already mentioned, the global transition prob- 
ability Y,(t¢) is obtained from the expression of Y, valid when the 
k,’s are independent of t , by performing the substitutions: 


ki = K; 5 | age gaat & B 
where it is convenient to define the new time variable T by the equa- 
tion 


r= | f(eydt. (15) 


We will identify the time variable t with the age of the microorgan- 
isms, each of the irradiated samples being assumed to be homogene- 
ous, so that all microorganisms of a sample are of the same age. The 
quantity 7 will be chosen as the age at the moment when the irradia- 
tion of the sample starts, + is in general different for each sample. 
The age is taken here not in an absolute sense, but with reference to 
the youngest. sample considered, 7 = 0 being taken as the moment in 
which the irradiation of that sample starts. 

The case which has heen treated up to now by various authors is 
the one in which the K;,’s are independent of 7. The meaning of this 
assumption has been discussed in Section 2. This is also the assump- 
tion that we consider here. We take k,.1 = 0, as before, that is, 
we assume that any change beyond the nth state does not affect 
the observation (the nth state is the final one). We put K; = KJ for 
4 =n, where K is a constant and J is the intensity of radiation, be- 
cause we assume that the transition intensity k; are proportional to 
J. Then it is known that (cf. e. g. Crowther, 1926): . 


Y,(t) =y(q,n) whereq=KJT, 
(16) 


r(a,n) = [‘[eeq/(n—1) Naa. 


¥(q, 7) is the incomplete gamma function. Using the symbols of K. 
Pearson’s tables (1934) we shall write: 


Y,(¢) =I(u,p) whereu=q/\V/n, p=n—1. 
Let Y(t) represent experimental values of the function Y,(¢). 
We consider the problem of determining the sensitivity 
k(t) =K Jf (t) 
from Y.(t). It will be seen that in order to solve the problem it is 
necessary to have several sets of values of Y,, each set correspond- 


ing to a different value of +, that is to a different age of the sample 
at the moment when the irradiation starts. The various sets of Y.(t) 
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must correspond to overlapping intervals of time t, that is, each age 
range must be covered by several samples. Otherwise the problem 
would have an infinite number of sojutions, since is unknown, and 
to each value of n corresponds one solution. In fact, the function 
Y.(t) is obtained from irradiation experiments, whereas the function 
y(q, ) is known, for instance through its values tabulated by K. 
Pearson. For each value of , the function q(t) is determined by the 
equation 

y(q,n) = Y.(¢), 
which has a unique solution gq = q(t). This solution can be obtained 
through the following numerical procedure; choose an arbitrary value 
q: of q, Pearson’s tables give the corresponding value y, of y for each 
assigned value of ; the experimental data give a value #, of t for 
which Y,.(t,) =y,. In this way for each value of g the corresponding 
value of t is obtained. A numerical or graphical differentiation of 
q(t) gives then the sensitivity according to the following relation (cf. 
equations 16 and 15): 

dq/dt = KJf (t). (17) 


The use of Pearson’s tables could be avoided by the following method, 
which is very complicated, however: Calculating the expression 


dlog Y’.(t) /dt, 


where Y’.(t) is taken as dY,(t)/dt and for Y,(t) the equations (15) 
and (16) are used, an integral equation of the type 


i ‘7 (t)dt=G(f, df/dt) (18) 


is obtained, where G is a rational fractional function of second de- 
gree in f and df/dt. A differentiation of equation (18) with respect 
to t gives a nonlinear differential equation of second order in f(f). 
The initial conditions for f(t) are complicated involving values of de- 
rivatives of f up to the mth order. However, for our purpose, the 
numerical method, as previously outlined, is completely sufficient. 
Let 7 = 0 and + = 7, be the two smallest values of + for which 
experiments are carried out. Consequently the range 0 < ¢ < 7, is 
covered by experiments on one sample only; therefore, the sensitiv- 
ity calculated in that range may be incorrect because of the arbitrary 
choice of the value of n. Before we show how this error may be cor- 
rected, we have to make some preliminary clarifications. We exclude 
the possibility of an infinite sensitivity at t = 0, because this would 
mean an instantaneous action of the radiation on the whole sample, 
which is in contradiction with experimental findings. Consequently, 
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since our range of time includes t = 0, we can expand f(t) in Mac- 
Laurin series of t, for sufficiently small values of ¢: 


f(t)=fot fit t---. (19) 


We assume that the experiments are pushed down to such a low value 
of +, that the series (19) converges for 7=7,. A particular meaning 
has the case f) = 0, because it requires that the sensitivity be zero 
at the beginning of the radiation (t = 0) and becomes apparent only 
as the radiation goes on. We have here a case in which the sensitivity 
to radiation is induced. by the radiation itself. There is, of course, 
still another interpretation of the condition f, = 0 possible, that is, 
that the microorganism exhibits a lack of sensitivity for some intrin- 
sic reason existing during a certain range of its age. Such a possi- 
bility can be cleared by an additional experiment which should start 
at an age slightly lower than 7,. We exclude these two possibilities 
from further consideration, so that we assume f, + 0. 

We show now how a possible erroneous determination of f(t) 
due to arbitrariness in the choice of n may be corrected. The func- 
tion f(t) represents the change of sensitivity. Let F(t) be its ex- 
pression determined from experimental values of Y.(t), by assuming 
a certain value N for the number n of effective events. Since F'(t) 
is not necessarily the correct expression of the change of sensitivity, 
the expansion of F(t) in power series has a general form 


F(t) =F @4+ Fywat 4+ --, (20) 
where 7 may be any real number, so that even the case F(0) = © is 
not excluded. Let n be the correct number of effective events and 
f(t) the corresponding correct expression of the change of sensitiv- 
ity, so that for f(t) the expansion (19) holds with f, #0. Since the 
pairs of quantities [f(t), n] and [F'(t), N] are derived from the same 
experimental data, they must satisfy the following relation (cf. equa- 
tions (15) and (16) ): 

7y(q,n) =7(Q,N), (21) 


where 
a=KI | f(tat; O=> KI [ r@at. (22) 


The derivative of equation (21) with respect to ¢ gives the relation: 
et gq" f(t) /(m—1) != e? Q™ F(t) /(N—1)! (23) 


Substituting here the expansions (19), (20) and taking into account 
the relations (22), we obtain for the lowest terms in equation (23): 
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[Fi/(N—1) !] [KJFi/ (i + 1) 3 Pha 
[fo/(m—1) 1] (KJf,) 1 t™ +... 
This relation is an identity, therefore: 
n=N(t+1), (24) 


which gives the correct number n of effective events on the basis of a 
fit of experimental data with any number N of events. This proce- 
dure will be illustrated in the following section. 


5. Application of the theory to experimental results. We apply 
now the theory to some experimental results on irradiation of Ascaris 
eggs with X-rays. G. Perthes (1904) was the first who observed an 
effect of X-rays on the process of cellular division in Ascaris eggs. N. 
M. Stevens (1909) carried out similar experiments with ultraviolet 
light and J. C. Mottram (1913) with beta and gamma rays of radium. 
The experiments have been repeated by many other authors since 
then (Holthusen, 1920, 1921; Schleip, 1923; Ruppert, 1924; Dognon, 
1926; Zuppinger, 1928; Liechti, 1929; Cook, 1939). The effect has 
been described as irregularities and prevention of development or 
change of rate of development. In counting the number of microor- 
ganisms which exhibit an effect of radiation, Holthusen (1920) dis- 
tinguishes two degrees of injury, one in which some form of worm is 
distinguishable, and the other one, in which the development reaches 
only a stage of an irregular cluster of cells. Other authors (e.g. Zup- 
pinger) do not differentiate between the degrees of injury in evaluat- 
ing the effect of radiation. From a theoretical viewpoint both proce- 
dures are admissible: the first involves two observable states and con- 
sequently two different values n, and 7m» of 7; in the second procedure 
one state only is observed, that one which corresponds to the smaller 
of the two values 7,, m,. For an application of the theory outlined 
in the previous sections the experimental results of A. Zuppinger 
seem to be more appropriate, because they cover various ages of the 
microorganisms and are among the most accurate that are available 
in this field. 

Figure 1 shows the results of 6 sets of experiments of A. Zuppin- 
ger, each one carried out with a different starting age. Three differ- 
ent intensities of radiation were used from 5.37 to 50 Roentgen/min. 
The experiment with 50 Roentgen/min. covering the range of ¢ from 
0 to 48 minutes was carried out in the earliest stage of development; 
this is the only experiment which has been discussed theoretically by 
other authors. By assuming a constant sensitivity k, A. Zuppinger 
(1928) obtained the closest fit with a number of effective events 
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n = 5; since the fit was not satisfactory he assumed a variability of n 
due to a non-homogeneity of the sample. The number ” was taken 
according to a Poisson distribution. The expression of the probabil- 
ity Y,(t) involves then Bessel functions and is much more compli- 
cated than in our theory. Nevertheless the fit could be improved only 
slightly by assuming the most probable value of n equal to 11. Re- 
cently N. Koyenuma (1943) obtained a very good fit for the above 
mentioned single experiment concerning the first 48 minutes after 
incubation. He did so by introducing an additional assumption, that 
the microorganism contains several sensitive volumes and that the 
injury becomes manifest when a certain minimum number of effec- 
tive events occurs in each of these volumes. His calculations involve 
series of Bessel functions and the fit is obtained with four sensitive’ 
volumes and n = 6 as the most probable number of events. Since the 
reaching of the two-cell stage requires 3 to 34 hours, whereas the 
whole experiment lasted 48 minutes it would have to be assumed that 
each cell contains all of the four sensitive volumes. Neither Zuppin- 
ger nor Koyenuma took into account the fact that the Ascaris eggs 
exhibit a strong variability of sensitivity during mitosis, as it has 
been positively established by Zuppinger himself as well as by C. 
Mottram (1913), H. Holthusen (1920) and A. Liechti (1929). 

The present analysis is not limited to the single experiment con- 
sidered in the calculations of A. Zuppinger and N. Koyenuma. It is 
based on several sets of data of A. Zuppinger on continuous irradia- 
tion of Ascaris eggs which have been conserved in aerobic conditions. 
Only a few sets of Zuppinger’s results were discarded; they concern 
some preliminary experiments, in which the irradiated samples were 
kept close to the X-ray tube, without protection against the heat 
emanated from the tube. [It was found by Zuppinger himself that 
the heat has a perturbatory effect. A. Dognon (1926) found also that 
beyond 22°C the sensitivity of Ascaris eggs to X-rays increases very 
rapidly with temperature. H. Holthusen (1920) emphasized also the 
same fact]. With this exception and with the exclusion of some data 
concerning special experimental conditions (Ascaris kept on ice and 
eggs in anaerobic conditions) all data of Zunppinger were used. 

The process has been assumed to be without recovery, in agree- 
ment with the experimental findings of A. Zuppinger and other au- 
thors. Zuppinger finds even that fractional irradiation produces more 
effect than a continuous one with the same amount of energy. Such 
increase of effect by fractional irradiation has been observed, how- 
ever, only over periods of days or weeks, and mostly on eggs which 
have been kept in anaerobiosis. A quite pronounced effect of recov- 
ery has been observed by E. V. Cook (1939), but again only over peri- 
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ods of weeks and exclusively at lower temperatures (5°C). According 
to Cook the capacity of recovery depends here on the very low metab- 
olic rate due to the low temperature, so that the X-rays produce their 
full effect long before the cleavage takes place. 

In a theoretical analysis of the results of A. Zuppinger it is there- 
fore justified to neglect the recovery and to consider the process as 
consisting of direct transitions only. The intensities of transition be- 
tween each two states will be given by expressions of the form of 
equation (17) as discussed in the previous section. As the first trial 
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Determination of the change of resistivity of Ascaris eggs to X-rays during 
mitosis. Experiments of A. Zuppinger (1928): sets A, D, &, F — intensity 50 
R/min., set B — 10.86 R/min., set C — 5.37 R/min. The experimental data 
are the same for both curves q/J and have common abscissae; for the curve 
q/J — 2 they are indicated with dots only. The interval between S and T repre- 
sents the time within which practically the full transition from the single to the 
two-cell stage is accomplished (90% single cells at S and 1% at T). For further 
details, reference is made to the paper by A. Zuppinger (set A — p. 663, Fig. 18 
right; sets B and C — p. 677, Fig. 19, 50 and 70.7 cm distance of the sample from 
the tube; sets D to F — p. 663, Fig. 13 left). The percentage of eggs injured by 
causes independent of radiation has been taken into account in the same way as by 
A. Zuppinger. The large dispersion in the sets B and C is probably an effect of 
the heat irradiated from the tube. The curve R — 2 represents the instantaneous 
resistivity in Roentgen. 
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value of the number of effective events we choose the one considered 
by A. Zuppinger: n = 5. From experimental data we calculate q/J 
as a function of the time t by means of Pearson’s tables of the incom- 
plete gamma function. These values are plotted in the diagram in 
minutes per Roentgen (curve g/J — 5). Each of the six experiments 
leads to a different curve. For convenience of drawing, these curves 
have been shifted upwards by various amounts, which are immaterial 
for our conclusions since we are interested only in the slopes of these 
curves. (These amounts are in R/min: set A — 0; set B — 0.13 and 
0.17; set-C — 0.35 and 0.74; set D — 0.38 and 0.94; set EH — 0.76 and 
1.7; set F — 1.1 and 2.08, where the first figures refer to the curve 
q/J — 2. and the second to the curve g/J — 5). The differentiation of 
q/J with respect to the time gives the function Kf(t) according to 
equation (17). From this the quantity n/[Kf(t)] is calculated. If 
the sensitivity were constant, n/(Kf) would represent the “average 
life’ of the microorganism in minutes under an irradiation of 1 Roent- 
gen per minute, that is, the resistivity of the microorganism im 
Roentgens. This is a consequence of equations (2) and (13). In the 
present case, however, f(t) is a function of time. Therefore we will 
call the quantity n/ (Kf) the instantaneous resistivity (curves R). It 
is seen from the diagram that assuming 5 effective events, the resis- 
tivity of the microorganism is zero at t — 0, that is, its sensitivity is 
infinite, which is impossible on the basis of the discussion of Section 
4, We use the method of that section to correct the number 7 and con- 
sequently to obtain the correct resistivity curve. Using the symbols 
of Section 4, we have: N = 5, and 7 is determined from the exponent 
of the lowest term in the expansion 


a7 I = {KF s/ (4+ 1)yt +... , (25) 


which is obtained from equations (20) and (22). The first term of 
expansion (25) is fitted to the curve (q/J — 5) between t = 0 and 
t= 20. Such a fit is easily done by using the elementary formula 


t+1=t,¢5/u, (26) 


where q, and q’, are the values of q/J and of its time derivative for 
t = t,, t, being in our case any value of ¢t between 0 and 20. These 
calculations give i = —0.6 . Consequently from equation (24) we have 
the correct number of effective events: n =2. Now, with n = 2 the 
calculations are carried out in the same way as before with n = 5. 
The results are shown in the diagram as the curves (q/9 — 2) and 
(k — 2). It should be pointed out, that the data available do not make 
it possible to carry out the calculations very accurately, and the re- 
sults obtained should be interpreted accordingly. The decrease of re- 
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sistivity R from 1,000 to about 250 Roentgen during the cellular divi- 
Sion is in agreement with the results of C. Mottram (1913), H. Hol- 
thusen (1921) and A. Liechti (1929) which estimated an increase of 
vulnerability of Ascaris eggs to X-rays from 2 to 8 times. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Theoretical studies have been carried out to determine the effect of 
the route of entry and the temporal sequence of dosage on the response. 
On the basis of continuity equations applied to a schematized system, 
equations relating the dosage and response! time have been derived for 
some special cases, and the form of this relationship is discussed. 


The effect of any material taken into the system of an organism 
is in general very complex and some of the effects produced mutually 
interact. We shall consider the relatively simple case where the re- 
sponse is not appreciably influenced by the other effects produced. 

We shall use the following abstraction. The material, whose ef- 
fect is to be observed, enters some part of the organism, is transported 
to the blood and subsequently is enabled to pass into a particular re- 
gion where the response is produced. We shall assume that the func- 
tion of one such region dominates the outcome. 

Let the amount of the chemical in the organ of entry at time t be 
M(t), in the blood—m(t), and in the sensitive organ—a(t). Let the 
rate of entry into the blood, regardless of the path of entry be R(t). 
If ym is the rate of transfer of the substance into the sensitive region, 
and bm is the rate of transfer of m into any other region plus any in- 
activation in the blood, then the value of m(t) is determined from 
the continuity equation 

dm 

7 R(t) —bm— ym. (1) 
We shall assume, as a first approximation, that b and y are constants 
thereby neglecting the possibility for a threshold of elimination as 
well as the reverse flow from the various organs. 


1 This paper is based in part on work done for the Office of Scientific Research 
and Development under Contract NDCre-132 with the University of Chicago. 
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Similarly, if fx is the inactivation rate in the sensitive organ, 
we have for the determination of x , the expression 


If the return flow from the various regions is not neglected, the 
above equations are replaced by a set of equations involving all the 
regions. Let the subscript 7 refer to the ith region, exclusive of the 
sensitive region, in which there is a concentration c; of the substance 
administered. Let a; be the partition coefficient for cellular fluid of 
a region with respect to blood multiplied by the blood volume. If 
8; is an inactivation rate per unit concentration and 0; is the product 
of the equivalent permeability and the effective area multiplied by the 
ratio of the volume of the region to that of the blood, then we have 
the equations 


dm 

—=- R(t) —> bi(m— aici) , (3) 
dt i 

de, b 

apne i (m— aic;) —— BC; , (4) 


together with an equation in x similar to equation (4). 

If the response to be studied is graded then we may suppose that 
the degree of response is some simple function of the quantity x. A 
more complex relationship will be considered subsequently. We shall 
consider an all or none type of response on the assumption that the 
quantity of the material in the sensitive region determines the re- 
sponse. Let «; be the threshold value of x for the response to be pro- 
duced in the jth individual of a population, and let be the median 
value of the #;’s. Let ¢ be the time required for the response to be 
observed after the time T, for which x(7) = #, together with any 
other constant time delays. It will be understood that t is properly 
accounted for, so that it need no longer be written explicitly. 

The individual thresholds, #;, are distributed in some manner 
over the population. Empirically it is found that for many responses 
the distribution is made very nearly a normal error function if the 
logarithm of the variable is used. Let p(z) be the normal error func- 
tion and define P(y) by 


P(y) = [p@ad. (5) 


Then, if 2 is the standard deviation, we have assumed that 
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1 
y =~ log x#/z, (6) 
4 


as & was defined as the median value of the @,’s. Since x is a func- 
tion of ¢, P is also a function of t as well as of some measure of the 
dosage. 

We shall now consider the effect of several modes of entry. 

(a) Ifthe rate of entry into the organ of entry is constant (e.g., 
constant injection, constant breathing of material coming in through 
the lungs), then the rate will be roughly proportional to the external 
concentration C. Thus R = H(C—m) =HCifC>>wm. If Cisa 
known function of time C(t), it is still possible to solve equations (1) 
and (2), with the initial conditions m =2—=0fort=0. 

For this case the response is observed at a time J + t,, where 
T is determined from «(T) = &, x(t) being determined from equa- 
tions (1) and (2) with R = HC. This would imply an exposure to 
C for a time 7, that is, t units of time less than that necessary to 
produce the response. But according to the above assumptions, the 
response has already been set in motion; thus, continuation of the 
exposure can, at most, increase the per cent which responds. In case 
the exposure is stopped at 7’, or t; seconds before the response is ob- 
served, then if HC is constant, the relationship between T and C is 
determined from 


[OEP 0 — 4) £-— (8 —b —»)yHC (1 — €8") 
aaa y BHC (e-nT — e-B7) , 


For small T and large C , CT? = 2%/yH,, and thus the quantity CT 
for median response becomes large for small T or large C. If the 
exposure continues to JT + t,, this effect is only accentuated. If, how- 
ever, exposure is only for a time 7, chosen so that after a total time 
T + t, a median response is obtained, then C7; is a constant for small 
+. The quantity Cr may remain approximately constant over a con- 
siderable range of values of z. Similarly for not too small or too large 
T , CT is approximately constant. 

For long times or small C, CT ~ Cr becomes large. If C = 
f(b + yv)/yH, then the response never occurs in more than half the 
population, and generally will not occur in any individual for small 
enough C. 

(b) Ifa dose D is impulsively given into some region and if 6M 
is the rate of entry into the blood stream, and eM is the inactivation 
rate in the organ of entry, then 

» aM 


ee M, 
=r (6 + e) (8) 


(7) 
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where M =D fort =0. From equation (8), we have 
<= De)", (9) 


and thus 
R(t) = 6M = 6De-“®)*, (10) 


In this case equations (1) and (2) are solved with m = x = 0 for 
t=0. Solving these equations and setting x =z fort=—T, then T is 
determined from 


(0 + yd —6) B= ODy( (11) 


e-(6+e)T — e-BT e- Oy) TF — “—) 
p—d—e Bib 

If D is too small, evidently no T satisfies equation (11), thus no re- 

sponse will be obtained for small enough D. 

If D is very large, T very small, then T = \/2%/éyD and thus T 
increases very slowly as D decreases. The above result cannot be 
readily verified because of the variability of f,. 

(c) If the material is injected directly into the blood stream at 
a rate R which may be a function R(t), then equations (1) and (2) 
can be solved directly with m= «—0,fort=0. 

The results, for constant R , are the same as for case (a), if HC 
is replaced by R. In this case again the result depends upon whether 
the injection is continued until the response occurs or continued until 
some time 7 , such that some time later, T + t, —7+, a median response 
is observed. In the first instance the total dose D = RT becomes in- 
finite for large R. In the later case the dose D = Rr remains finite. 
For small R , D becomes infinite and if R < 6(b + y)/y no median re- 
sponse occurs. To determine D(r), first solve for 2(7) as in equation 
(7), and for m(z) from equation (1), and determine x(t) for t > + 
from equations (1) and (2) with R=0 and m= m(r) and x = x(r) 
for t= 0 as initial conditions. Differentiating x(t) with respect to t 
and setting the result equal to zero, the time 7'(7,R) can be obtained. 
For this value of the time the quantity x has reached its median 
threshold so that «[7(7,R)] =. This is an implicit relationship be- 
tween R and r-, or since D = Rr, we have the relationship D(r). 
Moreover since T is already determined, D determines + and T. 

(d) If the total amount D of material is injected very rapidly 
into the blood stream, then equations (1) and (2) are solved with 
R= 0 but with m= D,x=0 for t=0. This leads to the equation 


(8—b—y) €= yD (eo? — ef), (12) 


for determining the relationship between T and D. If D is too small 
there will be no response. If D is sufficiently layge then T = @/yD. 
For both cases (b) and (c) the threshold dose D,, can be obtained by 
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setting the quantity in parentheses on the right-hand side equal to its 
maximum value and solving this equation for D = D,, . 

Case (d) can be considered a special case of (c) for which R is 
large and 7 small, but such that Rr = D, the total dose. 

(e) If the critical organ is directly perfused with the chemical 
at a concentration C or if the organism under observation is unicellu- 
lar, then «(t) is determined from equation (2) with ym replaced by 
yC where y is the permeability times the effective area divided by the 
volume. In this case, if C is constant, C(7) can be solved for T and 
written as 


P= —t=—“Ioe(1-—=). (13) 


In each of the above cases if x is identified with the response 
measured, then x(t) is determined for the various cases by replacing 
T by tand # by z. 

Cases (c), (d), and (e) are determined by the three parameters 
&, (6b + y) and &/y. With the additional parameter H , case (a) is 
also specified, whereas e and 6 given in addition to the first three 
parameters, case (b) is also determined. 

For each of the above cases, we can determine, by equations (5) 
and (6), the relationship P[y(D,t)]. This determines the per cent 
response within the time t if the dose is D. 

If t is fixed, then y of equation (6) is linear in the logarithm of 
D, the slope being 71, and the intercept depending upon the value of 
t. If D is fixed, then y can be considered to be a function of ¢ with 
D as a parameter. In the cases which can be approximated by equa- 
tion (14) then, for not too large t, y is approximately linear in the 
logarithm of t , the slope being again 4, and the intercept increasing 
with D. For larger t¢, the curves become concave downwards, ap- 
proaching a limit. On the other hand, if equation (14) cannot be used, 
then the slopes of y of log ¢ for not too large ¢, may no longer be /*. 
But in either case the results are complicated by the variability of t. 

In equation (12) if 6 and (b + y) are far from equal and if a 
is a number equal to the larger, then to a fair degree of approximation 
the time T is given by equation (13) with # replaced by a and C 
replaced by D. From this, T approaches infinity as D approaches 
D,» = a&%/y. Experimentally equation (13) is almost indistinguish- 
able from 

Ne Rea SE, 6 ae LS ar (14) 


if x = £/y and n= 1, especially if t, is slightly changed (cf. J. Ipsen, 
1941). Furthermore, even if f and (b + y) are comparable, equa- 


126 RESPONSE TIME AND DOSAGE OF A DRUG 


tions (12) and (14) are experimentally almost equivalent, the prin- 
ciple difference being that in equation (12), T cannot exceed some 
large value 7, for D = D, whereas in equations (13) and (14) 
T,, = ©. This, however, is experimentally difficult to determine. 

Equation (7) is similarly almost equivalent to (14) whereas 
equation (11) may not be. If any one of the quantities 6, (b + y); 
and (6 + «) is considerably larger than the others, then equation (11) 
may be approximated by equation (14) for m= 1, with a new value 
of t. which depends slightly upon D and which reduces very slowly 
to t, for very large D. For practical purposes, this means that the 
observed t, may include a time delay which is only roughly constant. 
If, for example, all three quantities are equal to a, then instead of 
expression (11), we have 


2% = byDT2e"? (15) 


In this case, with an adjustment on D,, , a fairly satisfactory ap- 

proximation is given by equation (14) with nm = 1/2 and «x = V/2z/y6 
except for D very near D,, . 
From the above results it would be expected that the relationship 
between dose and response time could be represented in many in- 
stances by equation (14) in which nm = 1 and in most instances if n 
is some fraction. However, this would be based on the assumption 
that the instantaneous value of « determines the response. We shall 
now consider a more realistic picture. 

Let the critical organ normally produce the total amount F of 
some substance which is the quantity to be measured or which is ade- 
quate to produce the response to be observed. The quantity F' will, 
in general, be determined from a continuity equation which has the 
form 

ab G F 

7 oF, (16) 
where G and g are, in general, functions of F’, as well as the concen- 
trations of many other substances. Even if G and g are assumed to 
be constant as a first approximation, any impressed change in F from 
its equilibrium value F, = G/g results in a spontaneous return to that 
equilibrium. Equation (16) is the simplest equation representing the 
property of homeostasis. In particular G or g must depend upon 2. 
If F is produced outside the critical organ, then G is the rate of entry 
extrapolated to F = 0 while gF is the return flow due to F + 0 plus 
any inactivation. In this case it would be expected that G may be 
fairly independent of x whereas g in a first approximation would be 
given by g = go + rz, that is, the removal of F due to x would take 
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place at a rate proportional to x and to F. 

If F is produced in the critical organ, then either G or g may 
depend upon x. If the effect of x is to reduce F then g(x) may ina 
first approximation be of the form g = g, + 7x over a wide range 
of x, whereas G(x) in the form G) — k’x would be expected to hold 
only for a limited range of values of «. If x acts to increase F then 
to a first approximation g would probably be independent of « , where- 
as in this case G could be represented by G = G, + kx as a first ap- 
proximation. 

If g is aconstant and G(x) = G, + ka then for case (d) we have 
an equation similar to equation (11). The equation for D(T) in this 
case is obtained from (11) if f is replaced by (b + y) and @ is re- 
placed by (F — F.)/k, F being the threshold for the response to be 
observed. 

If, on the other hand, we consider that x acts to decrease F', the 
result would be the same if F < F, and F, — F << F,. This latter 
inequality is rather unlikely. It is more likely that the function of 
the critical organ must be reduced to at least one-half or even to a much 
smaller function, so that the approximation G = G, — k’x would be 
of little value. Even if for each physiological unit G could be repre- 
sented by y = G, — k'x for x < G,/k’, y= 0 for x > G,/k’, the fact 
that not all would have exactly the same parameters results in F grad- 
ually approaching zero. In writing equation (2) for x we have as- 
sumed essentially an infinite internal diffusion coefficient. To a first 
approximation, removal of this restriction simply subtracts a constant 
from «. Thus F(x) for t ~ © could be represented by a curve hori- 
zontal at x = 0, having an inflection point, and which then gradually 
approaches zero, though being zero for some very large x. But if the in- 
terest is only in the response produced when F = F’, with F, several 
times F , then the shape of F(x) for near F, may be of no interest. 
Furthermore in a quantized response only the time at which F = F 
is of interest. This tends to reduce the effect of changes in the form 
of F(x) on the D(T) relationship. Nevertheless, it may be that the 
final relationship between D and T is no longer approximately given 
by equation (14) forn =1. 
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